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A statistical multistream description of quantum plasmas is formulated, using the Wigner-Poisson 
system as dynamical equations. A linear stability analysis of this system is carried out, and it is 
shown that a Landau-like damping of plane wave perturbations occurs due to the broadening of the 
background Wigner function that arises as a consequence of statistical variations of the wave function 
phase. The Landau-like damping is shown to suppress instabilities of the one- and two-stream type. 
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I. INTRODUCTION 

It has recently been pointed out [jj that the persistent 
trend towards increased miniaturization of electronic de- 
vices implies that quantum effects will become impor- 
tant also for certain transport processes, for which so 
far classical models have been sufficient. An example of 
such a generalized transport equation, in the form the 
Schrodinger-Poisson equation was analyzed in Ref. 0]. 
This analysis of a quantum plasma is based on the hy- 
drodynamic formulation of the Schrodinger-Poisson sys- 
tem, where macroscopic plasma quantities such as den- 
sity and average velocity are introduced. However, the 
analysis does not take into account statistical (or kinetic) 
effects associated with the finite width of the probability 
distribution function. Kinetic effects are well-known in 
plasma physics, where they may lead to the phenomenon 
of Landau damping. 

The possibilities of using a general approach based on 
the Wigner function 0, 0| was commented upon in Ref. 
0, but only a simpler approach based on macroscopic 
quantities was used. Obviously, in doing so the possi- 
bilities of Landau-damping like effects are lost. In fact, 
the possibility of obtaining Landau damping is also men- 
tioned in Ref. 0, although in connection with a possible 
generalization to the multi-stream case, in accordance 
with the classical picture of Dawson Particular at- 
tention was given to the classical one- and two stream 
instabilities in a cold plasma and it was shown that the 
main quantum effect on the wave propagation could be 
characterized as a generalized dispersion. 

However, recently much attention, within the nonlin- 
ear optics community, has been devoted to effects of par- 
tial wave incoherence e.g. in the form of phase noise on 
a constant amplitude wave 0E10]. In particular, it has 
been shown in Ref. Q, where the Wigner transform was 
introduced as a means to study the modulational insta- 
bility of an optical plane wave, that the phase noise gives 
rise to a Landau-like damping effect on the one stream 
modulational instability. 

It is the purpose of the present work to generalize the 
analysis made in Ref. by analyzing the properties 
of the one- and two-stream instabilities in a quantum 
plasma using the Wigner formalism and including the 



effect of phase noise developed in Ref. 0. The results 
clearly show the suppressing effect on the instabilities due 
to the Landau-like damping effect caused by the phase 
noise of the Wigner function. 



II. QUANTUM STATISTICAL DYNAMICS 

In non-relativistic many-body problems, the Wigner 
transformation is a useful means to derive equations de- 
scribing the quantum statistical dynamics of the system 
of interest. Thus, one is able to generalize the classical 
Vlasov equation to a quantum mechanical regime, in the 
sense that the dynamical equation for the Wigner func- 
tion describes particles moving in a self-consistent force 
field and in such a way that the evolution equation for the 
Wigner function takes the form of its classical analogue 
in the limit h — > 0. 

Haas et al. have considered the dynamics of a 
quantum plasma described by the nonlinear Schrodinger- 
Poisson system of equations: 
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where i — 1, TV numbers the electrons as described by 
pure states, with tpi being the wave function for each such 
state; (j)(x,t) is the electrostatic potential, while m and 
— e are the mass and charge of the electrons, respectively. 
The fixed ion background has the density no. Following 
Ref. 0], we have introduced the Klimontovich statistical 
average, denoting it by (■). The statistical averaging be- 
comes important when the wave function contains e.g. a 
stochastically varying phase 0. 

In Ref. 0|, the one-stream and two-stream models have 
been investigated and the dispersion relation for the two- 
stream instability was derived, showing an appearance of 
a new, purely quantum branch. We note that the anal- 
ysis presented in Ref. 0] is based on the hydrodynamic 
formulation of the system JQl , where macroscopic plasma 
quantities, such as density and average velocity, are intro- 
duced. However, this type of analysis does not take into 
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account statistical properties of the wave function that 
may lead to a broadening of the probability distribution 
function. In fact, such effects may give rise to a Landau- 
like damping both in the case of the single-stream and 
two-stream instabilities. 



In order to take the statistical effects into account, it is 
convenient to introduce the Wigner distribution function 
Wi(x, t;p), corresponding to the wave function ipi(x,t), 
as 



+ OC 



Wi{x,t;p) = — / dy exp(ipy/K)(ipf(x + y/2,t)ipi(x - y/2,t)) , 



(2) 
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which has the property 

f + OO 



dpWi{x,t;p) = (|Vi(M)r> 
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Using Eq. <|2j , Eq. Jlall can be formulated as a kinetic 
equation for the Wigner distribution, viz the Wigner- 
Moyal equation 



where the sine-operator is defined in terms of its Taylor 
expansion. Correspondingly, Eq. l|lbj) can be rewritten 

as 



dx 2 ~ e Q \ 2-t 



dp Wi - n 
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Clearly, an equilibrium solution of Eqs. Q an d CEt is 
4 = and Wi=Wio(p). 

In order to study the modulational stability of the sys- 
tem Q-© , we introduce a small perturbation according 
to 

Wi(x,t;p) = Wio(p) + Wi eicp[i(Kx - fit)] , (6a) 
<j>(x, t) = <j> exp[i(Kx - fit)} , (6b) 

where |Wj| <C \Wio\ and K and fl are the wave number 
and frequency of the perturbation, respectively. The fact 
that the background distribution is assumed to be 
only a function of p corresponds to the assumption of a 
plane wave function with constant amplitude, but with a 
stochastically varying phase, the characteristic properties 
of which are expressed by Wm(p). Linearizing Eqs. (@J 
and JU, we obtain 
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where (f> is the potential perturbation. Note that the fact 
that the unperturbed potential <f> is 0o = means that 
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Eliminating (j> in Eqs. we obtain the dispersion 
relation 
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= i [W M ( P + UK /2) - W i0 (p - KK/2)] (10) 
relation © can be written in the form 
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W M (p + HK/2) - W l0 (p - HK/2) 
p — Qm/K 



(11) 

Note that the pole p = Qm/K gives rise to both a prin- 
cipal part and an imaginary residue, as in the classical 
analysis of Landau damping in plasma physics. 

Let us now consider the cases of one-stream and two- 
stream plasmas. 



A. One-stream plasma 

The dispersion relation 11 111 reduces to 

+°° W (p + UK/2) - W (p - UK/2) 



1 = 
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p — Qm/K 
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where Wo = Wiq. For a one-component Wigner spec- 
trum with a deterministic phase, i.e. a monoenergetic 
beam, Wo(p) is given by 



W (p) = n S(p~p ) 



(13) 
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which corresponds to a monochromatic plane wave func- 
tion with constant amplitude and phase. Equation l(T2|l 
then yields 
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where vq = po/m and uj 2 = noe 2 /meo. The expression 
II 1511 is exactly the same as the one obtained in Ref . Q . It 
shows that quantum effects give rise to wave dispersion 
for short wave-lengths. 

Let us now assume that the phase <p(x) of the wave 
function ipo varies stochastically, and that the corre- 
sponding correlation function is given by 



l e -iip(x+y/2) e i<p(x-y/2)\ _ p -pT\y\ 



This corresponds to the Lorentzian spectrum 

Wo{p) = —-, rg- 2" > 

7T (P-Po) 2 +Pt 

and the dispersion relation Ill2t now yields 
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This result implies a completely new effect, a Landau- 
like damping due to the width of the spectral distribution 
describing the stochastic variation of the phase, i.e. due 
to the partial incoherence of the beam. Furthermore, the 
damping effect increases with increasing incoherence, i.e. 
with increasing pr- 



B. Two-stream plasma 

According to Eq. 11 111 , the dispersion relation becomes 
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For monochromatic beams with 

Wjo{p) = n j6(p-p j) ; j = 1,2 , 
we get from Eq. l(T9l) 
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where loz, = e nnj /enm and nni H-nn2 = ^o- If we follow 
Ref. [1| and consider the symmetric case where n i = 



"02 = n /2, poi 



~P02 



CI 4 - 1 + 2K 2 



- po, we obtain 

H 2 K 4 " 
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from Eq. iQTJ. Here we have introduced dimensionless 
variables according to 



CI = Cl/ujpo , K = pqK/c 



H = huj p om/pl . 



(23) 



Equation il22ll is identical to the result obtained from the 
hydrodynamical theory, as in Ref. Q. The solution of 
Eq. (St is 



CI 2 = ~ + K 2 



H 2 K 



! ± ^Vl + 8K 2 + AH 2 K e , (24) 



which implies O 2 < and concomitant instability if 

[H 2 K 2 -A){H 2 K 4 -AK 2 + A) < . (25) 
This condition can be written 



1 H 2 K 2 



(26) 



which reduces to the well-known two-stream instability 
result K 2 < 1 in the classical limit H — > 0. 

However, we infer from Eq. 12511 that the quantum ef- 
fect has a subtle influence on the instability. Equation 
ll25ll implies instability when the following condition is 
satisfied in (K, H) space, viz 



H 2 _(K) 



K 2 



1 
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K 2 



<H Z < — = H 2 + (K) . (27) 



A qualitative plot of this is given in Fig. (a similar 
figure and discussion was given in Ref. [J, but for later 
reference we present the figure and a discussion related 
to it). 

Figure ^ implies that when H = 0, instability occurs 
only for < K < 1. However, when H ^ 0, a more 
complicated picture emerges. In fact, as is seen from 
Fig. El the quantum effect plays both a stabilizing and a 
destabilizing role. For H > 1, instability occurs for all K 
such that < K < K + (H) = 2/H. Thus, for 1 < H < 2, 
the region of instability is increased, whereas for H > 2 
it is decreased as compared to the case H = 0. 

For H < 1, instability occurs in two -ftT-bands, viz 
< K < K W (H) and K {2) (H) < K < K+(H), 

( 1 2) 

where K_ : (H) are the two solutions of the equation 
1-1/K 2 = H 2 K 2 /A, i.e. 



K^(H) = JL (l + VT^JT 2 ) , 



(28a) 
(28b) 
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FIG. 1: Qualitative plot of the stability /instability regions for 
the two-stream quantum plasma, neglecting stochastic damp- 
ing. 



For all values of H < 1, this implies a larger range of 
unstable wave numbers as compared to the classical case 
H = 0. 

Let us now assume that the unperturbed Wigner dis- 
tributions have Lorentzian form, in analogy to the case 
of a one-stream plasma, i.e. 



n 0j 



PTj 



n (p-Poj) 2 +PTj 



J = 1,2 



(29) 



From Eq. ijlfljl we then obtain 



1 = 



[CI - (pox - ip T i)K/m] 2 - h K A / Am 2 

+ ^2 

[fi - (P02 - ip T 2)K/m} 2 - h 2 K 4 /4m 2 



(30) 



Following Ref. we consider the case when pox = 
— P02 = Po and not = no2 = no/2, while for the statis- 
tical broadening we assume pn = PT2 = Pt- Using the 
dimensionless variables given by (l23l) . we get 

(fl + iaK) 2 = l+K 2 + S^l±ly/i + 8K 2 + 4H 2 K e , 

2 4 2 (31) 

where we have introduced the relative broadening a = 
Pt/po- Thus, in the limit px — > 0, we regain the result of 
Eq. and Ref. 0]. However, in the previously unstable 
region we now obtain 



Im(fi) 



-aK 



i(l + 8/? 2 + 4ff 2 X 6 ) 1/2 
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H 2 K 4 



1/2 



(32) 



Again, the broadening a tends to suppress the growth, 
and the condition Im(O) > is now given by 



a < 



1 

K 



1 



1 + 8K 2 



4H 2 K 6 ) 1/2 



-2- KZ 



H 2 K 4 



1/2 



(33) 



In the classical limit H — ► 0, the region of unstable 
^-values is reduced to K < K c by the damping effect, 
where 



K,= 



1 



< 1 



(34) 



Clearly, for a > 1, no instability is possible for any K. 
Another illustration of this is the small- AT expansion of 
the growth rate, which reads 



Im(fi) ~ (l-a)K 



(35) 



The stabilizing influence of a in the general case of 
H ^ can be inferred as follows: 

Consider first the case of small K, while keeping 
H 2 K 2 /4 ~ ^(1), i.e. we investigate the growth rate close 
to the stability boundary. In this limit we obtain 



Im(O) 




H 2 K 2 



K 



(36) 



which clearly shows the stabilizing effect of the damping. 
In particular, the stability threshold is now given by 



(37) 



Qualitatively this implies a lowering of the upper thresh- 
old curve for small K and a concomitant decrease of the 
region of instability. 

Consider next the limit K 1, while still assuming 
H 2 K 2 /4 ~ ^(1), i.e. we examine the effects of the damp- 
ing on the narrow instability region, see Fig.Q] Introduce 
the notation 



Ah = 1 



H 2 K 2 



(38) 



The growth rate can then be written as 



Im(O) 



-aK 



I Ah [\-K 2 Ah] , (39) 



and the stability thresholds become determined by 



Ah 



1 



AK 2 



± 



1 



16K 4 



a* 



(40) 



When a = 0, we regain the previous limit curves Ah = 
and Ah = \/(2K 2 ). The effect of a nonzero a is to 
narrow the instability region and to terminate it at the 
finite wave number K = l/(2 v / a). For increasing a, the 
unstable region decreases and, as in the case of small K, 
we expect the instability to be essentially quenched for 
a > 1 . 
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III. DISCUSSION 

In this work, we have presented an analysis of a multi- 
stream quantum plasma, including the effect of phase 
noise in the beam wave function. As compared the fluid 
description of a quantum plasma used in Ref. pj, the 
present analysis is based on the quantum mechanical 
Wigner formalism. The phase noise, or partial incoher- 
ence, of the beam wave functions is shown to give rise 
to a Landau-like damping effect, which tends to sup- 
press the instabilities occurring in both the one- and two 
beam cases. The damping rate increases with increas- 
ing degree of incoherence as expressed by the width of 
the probability distribution function for the phase noise. 



The physical origin of this damping effect is the non- 
coherent properties of the beam wave function as op- 
posed to the wave-particle interaction characteristic of 
the conventional Landau damping. The new Landau-like 
effect is not a true wave damping, but a conservative 
rearrangement of the spectrum of the beam wave func- 
tion. This phenomenon has recently attracted consider- 
able interest, both theoretically[5, 6] and experimentally 
Q, within the area of nonlinear optics, where it has been 
shown to suppress the modulational and self-focusing in- 
stabilities (a&LLlI, e.g. for optical beams in nonlinear 
photo-refractive media. The present work is the first at- 
tempt to extend this theory to a quantum plasma. 
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